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ABSTRACT 


This study examines 48 Ml tank components for possible application of an opti¬ 
mum age replacement policy. The purpose is to support a broader study associated with 
the Reliability Centered - Inspect and Repair Only as Necessary (RC-IRON) program. 
The program provides depot level maintenance to tanks transferred or retrograded from 
Germany to the United States. An optimal age replacement policy reduces the number 
of failures while minimizing the cost associated w'ith failure by replacing some older 
components before they fail. The component data for this analysis was drawn from the 
Field Exercise Data Collection (FEDC) at the National Training Center (NTC), Fort 
Irwin, California. 

This thesis begins with a discussion of a methodology for determining an optimal 
replacement time. Distribution analysis is peiformed on component lifetimes as well as 
delay and repair times due to failure. The various costs associated with failure are esti¬ 
mated. The application of an age replacement policy was found to be beneficial for a 
few components and only when they had a high down-tirnc cost. A graphical procedure 
is used to show sensitivity of the optimum policy to changes in cost. Component simu¬ 
lations arc performed to pretest the results of a proposed maintenance policy. A six 
component system is simulated to demonstrate how the components could be tied to¬ 
gether for later system level analysis. Although this study deals with the Ml Abrams 
tank, the methodology and procedures detailed may be applied to other systems with 
components that wear out. 
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THESIS DISCLAIMER 


The reader is cautioned that all computer programs developed in this thesis research 
may not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of logic and computa¬ 
tional errors, they Cannot be considered validated. Any application of these programs 
without additional verification is at the risk of the user. 







I. INTRODUCTION 


A. PROJECT 

The Ml Abrams Main Battle Tank is a centerpiece of the modern battlefield and 
may not be replaced for quite sometime. One of its orginal design features was that the 
components, not the entire system would be repaired at depot. The Army has indi¬ 
cations that the burden of sustaining the tanks has grown as they have aged. Because 
of this growing maintenance the Army leadership would like to know if it is economical 
to use depots to identify a tank's condition and perform maintenance to extend service¬ 
ability. 

The Deputy Chief of StafT for Logistics (DSLOG) tasked Headquarters, U.S. Army 
Materiel Command (HQ, AMC) to determine an optimal point in the life of the tank 
when a Reliability Centered - Inspect and Repair Only as Necessary (RC-IROX) pro¬ 
gram, should be applied and if the program can be improved. The life of a tank is 
measured in miles. Determining the optimal point for acceptance into the program is 
based upon operational conditions and economic analysis. The actual depot inspection 
and repair procedure has been established based upon previous programs. The overall 
viability of the program could be enhanced if improvements to the current depot imple¬ 
mentation were found. HQ, AMC tasked Tank Automotive Command (TACOM) 
project management responsibility and U.S. Army Materiel Systems Analysis Activity 
(AMSAA) with analysis support. Descriptions are provided for the Ml tank in Ap¬ 
pendix A and maintenance levels in Appendix B. 
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B. BACKGROUND 

Analysis by TACOM and AMSAA will be conducted initially over the next two 
years. To enhance the analysis efTort a basic hardware validation test was established. 
This test consists of transf'ering or retrograding 60 Improved Ml (IPM1) tanks from the 
1st Infantry Division (Forward) in Germany to the National Training Center (NTC), 
Fort Irwin, California. At random, 14 of these tanks were selected and diverted to 
Anniston Army Depot, Anniston, Alabama, for RC-IRON, programmed to cost S95,000 
per tank, before being delivered to the NTC as the sample group. The remaining 46 
tanks were only subjected to standard deprocessing treatment. These tanks became the 
control group. The tanks arrived at the NTC in the first two months of 1990 and started 
to be used in exercises during the summer. The NTC was selected because the most 
mileage can be accumulated there and a data collection program is already in progress. 

C. OBJECTIVE 

The purpose of this study is to aid the RC-IRON project by focusing on the optimal 
age replacement policy for single components. A system level analysis is also breifiy 
considered. Them re 4S candidate components that AMSAA indicated were of prime 
interest. TACOM and AMSAA selected these parts because they are costly due to the 
high value of their cost per unit and the frequency of repair. In this thesis we will at¬ 
tempt to find an optimal maintenance policy for components among the suggested 48 
for which such a policy makes sence. The methodology is defined in Chapter II. Data 
used to examine the hS components is described in Chapter III. An initial screening of 
components is conducted in Chapter IV. Costs associated with this study are detailed 
in Chapter V. The components for which age replacement policies are appropriate, will 
be determined in Chapter IV. Based on the data in Chapter III and the costs from 
Chapter IV, optimum replacement mileages are estimated for the appropriate compo- 
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nents in Chapter VI. Insights into system level analysis are also included in Chapter 
VI. 

D. MAINTENANCE POLICIES 

The ability of an armor unit to perform its mission is dependent on its tanks. A key 
factor is the quality of operation and availability which are primarily facilitated by 
maintenance. The two major types of maintenance are preventive and corrective. Cor¬ 
rective maintenance is performed by repair, replacement, or overhaul of equipment after 
it has failed. Often preventive maintenance is applied to extend service life or reduce the 
probability of failure (Ref. 1: p. 17). Due to the possible economic and operational 
benefits, this study concentrates on planned preventive maintenance to reduce the 
probability of failure. 

In particular, the preventative maintenance policy selected for study is the policy 
based on age (age replacement). This policy is implemented by making replacements 
either at the time of failure or after cp units of mileage. This is not to be confused with 
block replacement, where the policy is instituded by replacing a set of components in the 
tank at prescribed mileages k<p (k - 1,2,...) independent of the history of failures in the 
tank system. The advantage of block replacement is that it is easier to implement due 
to a decreased administrative burden. Management of the policy is simplified when the 
incident mileage need not be recorded. However, components are replaced more fre¬ 
quently than needed, under a block replacement policy, thereby leading to increased cost 
[Ref. 2: p. 158]. Thus, this study will focus on age replacement policies. The optimal 
age replacement policy minimizes long run expected cost per mile with replacement at 
a certain mileage (p*. The methodology for determining this cp * mileage is detailed in 
the next chapter. 
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II. METHODOLOGY FOR DETERMINING REPLACEMENT TIME 


A. REQUIREMENTS 

In order for preventative maintenance, under an age replacement policy, to be lu¬ 
crative, the cost of an unscheduled replacement at failure must be higher than the cost 
of scheduled replacement. A description and estimate of these costs may be found in 
Chapter V with C, representing unscheduled and C 2 scheduled maintenance cost. It is 
also necessary that the component life distribution have a failure rate that increases with 
mileage [Ref. 3: p.46]. It would not make sense to replace an item that does not age 
or that is improving with age. To guard against choosing a replacement policy that ac¬ 
tually increases costs by making replacements too frequently, the optimal maintenance 
policy is selected by minimizing the expected cost of repair per mile [Ref. 1: pp. 19-24]. 

B. DISTRIBUTIONS APPROPRIATE FOR EARLY REPLACEMENT 

If a component is not expected to wear out it would be ridiculous to replace it before 
it fails. In other words, if it is improving with age, or staying the same, leave it alone. 
A class of distributions which captures a particular notion of aging is the Increasing 
Failure Rates (IFR) class of distributions i.e., those distributions with increasing failure 
rate [Ref. 2: p. 159]. Both Gamma and Weibull distributions are IFR when their shape 
parameter a is greater than one. The Weibull distribution is widely used for reliability 
analysis and is expected to be the most appropriate for this analysis. Gamma and other 
distribution can be examined in later analysis. 

When the underlying lifetime distribution is a member of the two parameter Weibuil 
family with shape parameter a and scale parameter A , the density is given by 
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t > 0, 


( 2 . 1 ) 


M “ “ 2 ( 2 ,)-' 



with failure rate 


K0 = 0iX , / 2r 0. (2.2) 

When a > 1.0, the failure rate in Equation (2.2) is strictly increasing to infinity. As we 
will see, this property guarantees that a unique and finite optimal replacement age <p* 
exists. The larger « is the more wear a component exhibits over time. Thus, a for greater 
than 1, the larger it is, the more appropriate it is for the component to be included in the 
application of our maintenance policy. Excluded from consideration is the exponential 
distribution, a - 1.0, and Weibull distributions with decreasing failure rate (a < 1.0). 
To give the reader a feel for the Weibull distributions used in this study, Figures 1 and 
2 show densities and failure rates with different shape parameters. The distribution se¬ 
lected for illustration include, a - 1.2, 1.4, 1.6, 1.8, 2.0. For comparison, the scale pa¬ 
rameter ). is adjusted so that the expected lifetime is 2.0. In our analysis, lifetime is 
measured in miles, but we will use "time" and mileage interchangeably. 
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Figure 1. The Weibul! Density Function f(t) with £(A’) = 2.0 
Source: Uyar, O. [Ref. 11] 
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Figure 2. The Failure Rate of The WeibuII Distribution with £(A’,) — 2.0 
Source: Uyar, O. [Ref. 11] 


C. PROCEDURE 

A little background is now provided on the mechanics of finding an optimal re¬ 
placement age. Under an age replacement policy, the time of planned replacement is 
specified as <p* so that components are replaced at cp* if they have not already failed. 
If the sequence of component lifetimes can be modeled as independent and identically 
distributed (iid) with distribution F, then the times between replacement form a renewal 
process. Thus, under such a policy with replacement at t, the long run expected cost per 
unit time C(t) can oe determined from the following equations jRef. 3: p. 87]. 
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CM = 


(2.3) 


C, X Fjt) + C 2 x Fjt) 

[ F(x) dx 
*'o 


Where X has distribution F and F - I - F is the survival function. Note that C(t), 
in Equation (2.3) is the ratio of the expected cost of repairing one component and the 
expected time (mileage) between repair. 

The Weibull survival function is given by 

f<r w \ t> o, 

Fit) = (2-4) 

ll, t< 0. 


Inserting the survival function of Equation (2.4) above into Equation (2.3) leads to 
the long run cost function below. 


CM = 



- e-<*‘>') + C 2 e- lA, r 

[e-^^dx 

*'o 


(2.5) 


See Figure 3 r, r the cost function plotted using the five Weibull distributions de¬ 
picted in Figure 1, with cost C, = 5.0 and C 2 - 1.0. It can be readily seen, especially at 
the higher shape parameters a , that C(t) indeed has a global minimum referred to as 
optimal age replacement time q>*. A proof that cp* exists and is unique when the failure 
rate increases to infinity, as it does for Weibuil distributions a >1.0, is given in 
[Ref. 2: pp. 161-168], In Chapter VI, this formulation will be applied to component 
distributions to estimate <p* to see at which mileage point, in the components life, it 
should be replaced. The component failure distributions will be estimated from the data 
described in the next chapter. 
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Figure* 3. The Long Run Expected Average Cost Curves Avitli ETA') = 2.0 
Source: Uyar, O. [Ref. 11] 
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IV. SCREENING OF CANDIDATE COMPONENTS 


A. STATISTICAL AND OPERATIONAL SIGNIFICANCE 

Prior to any further analysis, components that had fewer than 20 failures were ex¬ 
cluded. This eliminated those components in which fewer than 1/3 of the vehicles ex¬ 
perienced a first time failure. They were excluded because the estimated probability 
distributions would be very suspect based on such small data sets. 

One of the candidate components was the power pack that is actually made up of 
four modules. It was necessary to view each module seperately because of the different 
characteristics they have. Only the engine module of the four met the 20 failure mini¬ 
mum. 

The tank track and road wheels were also eliminated from further consideration. 
These two items were identified in previous testing and use as having unacceptable wear. 
A contract was let in 1988 to produce a new track with a 300 percent increase in ex¬ 
pected life. The following table is a list of those components that were not eliminated 
from the original 4S components by the above screening. 
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NOMENCLATURE 


ENGINE STARTER, GAS 


NAME 


STARTER 


TRANSMISSION ASSEMBLY TRANSMISSION 


GRIP ASSEMBLY, CONTR 
GUNNER'S 


NOZZLE ASSEMBLY, LULL 


DISTRIBUTION BOX 


LINK ADJUSTING TRACK, 
RIGHT 


TURRET NETWORKS BOX 


ELECTRO-MECH FUEL 


PUMP, FUEL ELECTRICAL 


ELECTRONIC CONTROL AS¬ 
SEMBLY 


HUB, WHEEL ASSEMBLY 


SPROKET WHEEL 


PUMP UNIT, ROTARY 


POWER CONTROL UNIT 



GRIP 


NOZZLE 


DISTBOX 


LINK 


TNBOX 


EMFUEL 


EPUMP 


ECASMBLY 


HUB 


SPROKET 


RPUMP 


SIGHT, GUNNER'S PRIMARY 


THERMAL RECIVER UNIT THERMALREC 



IMAGE CONTROL UNIT 



LASER RANGE FINDER 


TURBINE ENGINE 


LASERRF 


ENGINE 


NSN 


2990-01-094-1377 


2990-01-136-1206 


2520-01-157-3745 


2520-01-202-9S65 


1015-01-076-6865 


1015-01-076-6739 


2910-01-124-9325 


2910-01-214-2640 


6110-01-169-5164 


2530-01-164-5805 


1015-01-076-6688 


2910-01 


2910-01 


2910-01 


2910-01 


2590-01 


2530-01 


3020-01 


4320-01 


1240-01 


1240-01 


1240-01 


1240-01 


1240-01 


1240-01 


075-4926 


080-9132 


083-3153 


232-9687 


154-6656 


063-5666 


065-6209 


073-4S29 


204-5765 


074-8969 


162-0367 


132-1693 


152-5344 


074-8947 


1240-01-246-1872 


1240-01-074-S946 


1240-01-149-8302 


2835-01-120-3674 


2835-01-216-8639 


Tabic 1. REMAINING COMPONENTS 
















































































B. RELIABILITY ANALYSIS OF THE COMPONENTS 


In Chapter III, recall that all the failures were recorded as interval censored or right 

I 

censored data. For each of the components, failure distributions were fit 
nonparametrically and parametrically. See Figure 3 for an example. The sample 
nonparametric cumulative distribution function is a step function which is calculated 
using Turnbull's nonparametric maximum likelihood estimator based on right and in¬ 
terval censored data [Ref. 4 pp. 169-173]. Such a procedure distributes probability 
among the censoring intervals and to the right of the largest censoring interval when the 
largest observation is right censored. The data was also fit parametrically to a Weibull 
distribution using the method of maximum likelihood. The fits were generally quite 
good, see Figure 5 supporting the Figure 4 example. The outliers that have low mileage 
and relatively high percentiles may be explained under the phenomina of infant mortal¬ 
ity. All the significant and appropriate component distributional fits and percentile 
plots, along with a table of parameters estimates and standard errors are detailed in 
Appendix C. A summary of these that have increasing failure rate indicated by the es¬ 
timated shape parameter a > 1.0 is given in Table 2. 


15 









CUMULATIVE PROBABILITY 
0-2 D.4 D.6 0.8 


STARTER 

WEIBULL CUMULATIVE DISTRIBUTION FUNCTION, N=113 



Figure 4. Starter Distribution Fit 
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WEIBULL PROBABILITY PLOT. N=113 



100 1000 10000 
MILES 


Figure 5. Starter Percentile Plot 


PARI' 

SI 1 API: a 

SCALE /> = W- 

NUMBER OF 
FAILURES 

START I R 

1.1999 

4<>S9.76. 

57 

TRANSMISSION 

1.4455 

7479.5 

29 

GRIP 

1.5539 

4309.9 

54 

NO//IP. 

! .3384 

4438.72. 

56 

DISTBOX 

1.5315 

6858.6 

32 

I INK 

1.2SO! 

10239 0 

.. 

22 


Table 2. DISTRIBUTION ANALYSIS RI SlT IS 
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V. COST ESTIMATION 
A. MAINTENANCE COST 

In analysis of this type, cost is a major factor in the evaluation process. To compute 
an optimal replacement interval the unscheduled and scheduled maintenance cost must 
be estimated. Generally the maintenance cost (COM) may be calculated using the 
breakdown provided in [Ref. 5: p. 3S3], 

COM = (C omm + omx "b Qmj "h Qmt *b Comp "b C om y+ C om p ) (5-1) 

with 

C omm = Maintenance Personnel and Support Cost 

C cm , = Cost of Repair Parts 

C oms = Test and Support Equipntent Cost 

C om , = Transportation and Handling Cost 

C )m/ = Cost of Maintenance Facilities 

C cn<t = Cost of Technical Data 

For this maintenance policy study, the cost of maintenance will be calculated based 
upon maintenance personnel and parts cost only. The other costs are assumed to be 
either negligible, compared to other uncertainties, or not relevant, or possibly sunk, for 
this analysis. The maintenance facilities and test equipment have already been pur¬ 
chased and are considered sunk cost. The transportation and data collection costs are 
difficult to ascertain at this time and should play a minor role in a component replace¬ 
ment policy. This will not be true in the system level analysis, since these costs, espe¬ 
cially the transportation cost, will play a significant role. The labor rates were computed 
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to be $104/hr in the field and S!75/hr at depot by the author [Ref. 6: p. 23] and myself. 
These figures were crudely calculated to obtain a feel for the labor cost and may be too 
high. We used figures from the base line cost estimate, maintenance allocation chart, 
and RC-IRON estimate cost to calculate these man-hour costs. In this study, it is as¬ 
sumed that the labor rate is S50/hr for both locations. This is a standard labor rate for 
many civilian repair shops in areas of the United States. The parts cost Table 3 were 
obtained from the current Army Master Data File. 


PART 

COST (DOLLARS) 

STARTER 

S 794.00 

TRANSMISSION 

S 139,998.00 

GRIP 

S 1,955.00 

NOZZLE 

S 944.00 

DIST BOX 

S 12,021.00 

LINK 

S 488.00 


Table 3. PART COST 


B. PENALITY COST TOR DOWN TIME 

In private industry the cost of down-time is found by estimating the cost of lost re¬ 
venue. The military docs not have a profit motive to fall back on. In this study two 
different levels of penalties will be developed for management consideration. They are 
the stand-by and float penalties, named after two possible Army actions. '1 he first pen¬ 
alty, float, is named after the Operational Readiness Float (ORF) which is designed to 
improve the readiness of combat units. Extra combat systems, float tanks, are kept at 
an intermediate support maintenance unit for exchange with a customer whose tank 
cannot be repaired in a specified time, The second and larger penalty cost is stand-by. 
This penalty is for a tank which is standing by and ready tc go in the event one of a 
unit's tanks fail prior to going on a critical mission. The penalties are recorded in dollars 
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and are based upon system costs. This will enable us to use them in conjunction with 
actual maintenance cost to determine the unscheduled and scheduled maintenance cost. 
The acquisition cost of the tank varies from year to year, for this study S3,000,000 will 
be used. This was based upon the cost for MlAl's from [Ref. 7: p. 32] Table 4, prin¬ 
cipally the 1990 figure. 


YEAR 

COST PER TANK 

1979 

S 3,390,000.00 

1984 

S 2,047,000.00 

1990 

S 2,977,000.00 

1991 

S 3,552,000.00 


Table 4. AVERAGE COST PER M1A1 1979-1991 


In the Abrams Base Line Cost Estimate [Ref. 8] the average annual sustainment 
cost tank for the IPM1 is S562.500 and M1A1 is S514,900. For analysis purposes we 
shall use S500,000. The Army has issued a life cycle estimate of 20 years for planning 
purposes. Using the 20 year life cycle and average annual sustainment cost the 
sustainment cost f the Ml will be S10,000,000 over its lifetime. 

The lower of tnc two penalties is the float. It is based upon the ORF action and 
calculated system aquisition cost, neglecting the sustainment cost of the float vehicle. 
The actual ORF cost would be higher because some sustainment cost would be incurred. 


Downtime Cost Per Day 
Downtime Cost 
Float Penalty Cost Day 


Down Time in Days 
*■ 365 Days ^ *' 


Aquisition Cost 


20 }ear Life Cycle Cost 

S 41 1 

(Downtime in Days) ( ———— ) 


) 


S 41 i. 

Day 


Day 


(5.2) 


Stand-by is the larger penalty. The tank it represents is standing by and cost are 
calculated proportionate to aquisition and sustainment cost. The annual cost of the 


















crew which would realistically also have to be standing by is not included 
[Ref. 6: p. 22]. 


Downtime Cost Per Day 
Downtime Cost 
Stand — by Penally Cost Day 


- ( 


Down Time in Days 
365 Days 

(Downtime in Days) ( 
S 1,781. 

Day 


x Aquisition + Sustainment Cost 
( on v„— r ~ ) 

S 1,781 


Day 


20 Year Life Cycle Cost 
") (5.3) 


C. REPLACEMENT COST 

If under an age replacement policy, a component is to be replaced before it fails, the 
cost of failure must be higher than the cost of scheduled replacement. The cost of failure 
can be in the form of cost, danger, or lost time. These costs will be refered to as C, for 
all unscheduled and C 2 for scheduled maintenance. The costs are calculated with the 
following linear relationships. 


C, - a(MTD) + b 
C 2 = a (MTTR) + b 


(5.4) 


with 

a - Penalty Cost Per Day 

MTD = Mean Downtime due to Delay in Days 

MTTR - Mean Time to Repair m Days 

b = Part Cost + Labor 

Labor = ( - r~r~ ) x (MTTR) 

Manhour 

The values used to compute C, and C 2 in Equation (5.4) are given in Table 5. 
MTTR figures in man hours weie provided by |Rcf. 6: pp. J2-3]. It is assumed that 
although the MTTR times are often for two mechanics, these times arc representative 
of the delay for scheduled maintenance. MTTR was converted from hours to days for 
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standardization. The labor cost in representative of the $50 per hour labor rate times 
MTTR. The labor rate and parts cost were detailed in Section B of this chapter. The 


replacement cost calculated with both pentalty type are in Table 6. The table also in- 

Q 


eludes the cost ratio — 


needed in the next chapter for the graphical replace- 


(C, - Q) 

ment interval section. The MDT for each component in days was obtained from fitted 
mean calculated in Section D and represents the unscheduled delay. 


PART 

MTD 

(DAYS) 

MTTR 

(MNHRS) 

MTTR 

(DAYS) 

LABOR (S) 

b (S) 

STARTER 

.76 

2.0 

.08 

S 100. 

S S94. 

TRANS¬ 

MISSION 

1.34 

6.6 

.28 

S 330. 

S 140,328. 

GRIP 

. .81 

1.3 

.05 

S 65. 

S 2,020. 

NOZZLE 

.63 

7.8 

.33 

S 390. 

S 1,334. 

DISTBOX 

.87 

1.8 

.08 

S 90. 

S 12.111. 

LINK 

1.40 

2.8 

.12 

S 140. 

S 628. 


Table 5. COST EQUATION INPUTS 


PART 

c, 

Q 

COST RATIO 

STARTER 

S 2.24S. 

S 1,036. 

- .85 

TRANSMISSION 

S 142.715. 

S 140,82o. 

-74.55 

GRIP 

S 3,462. 

S 2,109. 

-1.56 

NOZZLE 

S 2,456. 

S 1,922. 

-3.60 

D1 S'l BOX 

S 13,660. 

S 12.253. 

-8.71 

LINK 

S 3,121 

S 842. 

-.37 


Table 6. REPLACEMENT COST WITH STAND-BY PENALTY 
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PART 

c, 

Q 

COST RATIO 

STARTER 

S 1,206. 

S 927. 

- 3.32 

TRANSMISSION 

S 140,879. 

S 140,443. 

-322.12 

GRIP 

S 2,353. 

S 2,041. 

-6.54 

NOZZLE 

S 1,593. 

S 1,470. 

-11.95 

DISTBOX 

S 12.469. 

S 12,141. 

-37.37 

LINK 

S 1,203. 

S 677. 

-1.29 


Table 7. REPLACEMENT COST WITH FLOAT PENALTY 

D. DELAY AND REPAIR TIME DISTRIBUTIONS 

Distributions were Fit to delay and repair times. Results from this analysis were used 
in Section C, to compute maintenance cost for some components. These costs in con¬ 
junction with the estimated component failure distribution are used in estimating opti¬ 
mal replacement mileage in Chapter VI. This analysis may indicate that additional 
components should be eliminated from further policy consideration. The results of this 
section will also be used in the simulations in Chapter VI. 

The six delay and two repair data sets for the remaining candidate component were 
fitted to the lognormal distribution. The lognormal distribution was chosen because it 
seemed to model delay and repair times. Most but not all fits were good. It was decided 
to stay with this model for these time distributions because of the way the data was 
collected. It is human nature to use rounded time increments, such as a fraction of a 
day. In the following example (starter), the bulk of observations are at a half a day. 
Figures 6, 7, and 8 are examples of the histogram, 

cumulative probability plot, and percentive plot for delay times. The three plots and 
analysis table information for the six components are contained in Appendix D. 
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RELATIVE FREQUENCY 
D.2 D.4 D.6 


Figure 


STARTER DELAY 

LOGNORMAL DENSITY FUNCTION, N=57 



6. Histogram for Starter Delay 
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Figure 7. Cumulative Distribution for Starter Delay 
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STARTER DELAY 

LOGNORMAL PROBABILITY PLOT, N=57 



DAYS 


Figure 8. Percentile Plot for Sarter Delay 

The repair fits were made for two components for latter simulation use. These fits 
with their graphical counter parts to the above Figures are located in Appendix E. 
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VI. MAINTENANCE POLICY ANALYSIS 


A. COMPONENT OPTIMAL REPLACEMENT MILEAGE 

The component failure distributions of Chapter IV and cost results of Chapter V are 
now applied to Equation (2.3). This is accomplished by using the APL program in Ap¬ 
pendix F. The cost function C(t) was estimated by using simulation. This was accom¬ 
plished by using 8,000 pseudo random lifetimes, generated from a Weibull distribution 
with parameters estimated from the component data described in Chapter III. The 
simulated C(t) is within S.01 of the actual C(t). As an example, the cost function C(t), 
based on 1,000 pseudo random numbers, is plotted (Figure 9). The optimum replace¬ 
ment milage and coresponding minimum cost per mile were found for each component 
by minimizing the simulated cost function. These results are located in Table 8 and 9. 
Under the smaller float penalty no components are recommended for age replacement. 
If the higher stand-by penalty is adopted, the only components that are recommended 
for early replacement are the link, starter, and grip. The starter and grip are marginally 
recommended, because their replacement points are near the end of their useful lifes. 
It should also be noted that the replacement mileage should be rounded up as long as 
the optimal cost is not changed significantly. This is due to the very large cost of re¬ 
placing the component too early versus the relatively smaller increases in cost if it is re¬ 
placed to late. 
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PART 

COST PER 100 
MILES 

REPLACEMENT 

MILEAGE 

MAXIMUM LIFE 
(MILES) 

STARTER 

S 1.36 

19,668. 

20.000. 

TRANSMISSION 

S 82.92 

25,000. 

25.000. 

GRIP 

S 5.05 

11,972. 

12,000. 

NOZZLE 

S 2.45 

15,974. 

16,000. 

DISTBOX 

S 10.13 

20,000. 

20,000. 

LINK 

S 4.29 

29,921. 

30.000. 


Table 8. OPTIMUM REPLACEMENT COST & MILEAGE WITH FLOAT 
PENALTY 


PART 

COST PER 100 
MILES 

REPLACEMENT 

MILEAGE 

MAXIMUM LIFE 
(MILES) 

STARTER 

S 2.56 

18,865. 

20,000. 

TRANSMISSION 

S 84.61 

25,000. 

25,000. 

GRIP 

S 7.31 

10.000. 

12,000. 

NOZZLE 

S 3.75 

15.842. 

16,000. 

DISTBOX 

S 11.17 

19,985. 

20.000. 

LINK 

S 10.07 

16,384. 

30,000. 


Table 9. OPTIMUM REPLACEMENT COST & MILEAGE WITH STAND-BY 
PENALTY 


B. GRAPHICAL DETERMINATION OF REPLACEMENT INTERVAL 

The determination of the optimal replacement mileage can also be made graphically, 
'lire method which uses a scaled total time on test plot is described in 
[Ref. 1: pp. II 3-116) by Barlow. The main advantage to this procedure is that sensi¬ 
tivity to changes in cost is readily apparent. Another important feature is that other cost 
ratios may be analyzed very quickly. We will use this method on the components re¬ 
commended for the application of a replacement policy in the previous section. For 
these components a scaled time on test plot is given in Figures 11, 12, and 13. The costs 
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ratio are obtained in Chapter V. The procedure for obtaining the scaled time on test 

plot may be found in Appendix G along with an APL program for assistance. To read 

C } 

the plot, a line is drawn from the cost ratio — —- —-r on the horizontal axis to the 

V^l ^2) 

tangency point on the time on test curve. From the tangency point one may vertically 
drop down and read the value of the cumulative failure distribution evaluated at the 
optimum replacement milage. The optimal replacement mileage may then be read from 
the plot of the cumulative distribution in Chapter VI and Appendix C. 

If the scaled time on test function is rather flat in the vicinity of the tangency point, 
then the optimum is not sensitive to small changes in the cost ratio. If a new cost ratio 
is of interest, then a new line may be drawn which is tangent 10 the total time on test 
function to find a new optimum. Now we have the convenience of not having to per¬ 
form more calculations or rerunning programs for different costs. 








SCALED TIME ON TEST SCALED TIME ON TEST 

D.2 0.4 as O.B 1.0 .0 0.2 0.4 0.6 O.B 


SCALED TIME ON TEST PLOTS 



COST*RATIO >< CUMULATIVE PROBABILITY ' 


FOR THE GRIP 



-i -i 


Figure 10. Starter & Grip Scaled Time on Test Plots 
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Figure 11. Link Scaled Time on Test Plot 

C. COMPONEM ; SIMULATION 

The operation of the link and starter will now be simulated individually to pretest 
the results of the proposed maintenance policies. Previously gathered information was 
used as inputs to the simulation. Since the main data set in this analysis was from the 
NTC and that is where the tanks are going to, we will use the NTC operating tempo for 
accumulated mileage. That is, the functioning components will accumulate mileage at 
rate equivalent to an average NT C tank with no other component induced down time. 
The actual program was coded in SIMSCRIPT 11.5 using discrete-event methodology 
represented by the process in Appendix II. The simulations were run for two year and 
twenty year intervals for 5,000 repetitions a piece. In the simulation, tanks are 
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retrograted after two years and the life cycle of a tank is taken to be twenty years. Re¬ 


sulting failures, cost, and availability for 20 year runs are located in Table 8 along with 
their deviations. No parts were replaced for either component in the two year runs. The 
starter had 156 parts preventatively replaced or approximately 3% of the tanks had the 
policy applied during their lifetime. A larger number occured for the link 4270 or 85% 
experienced preventative replacements during a life cycle. 

The performance of the maintenance policies could be improved some, by not al¬ 
lowing replacements to occur towards the very end of the 20 year life cycle. This would 
serve to reduce some of the maintenance cost where the full benefits would not be real¬ 
ized. This simulation is a specialized version, with more component statistics taken, of 
the program discussed in the next section. The comments for either program apply to 
the other. Verification will be discussed in the next section. 



WITH OUT POLICY 

WITH POLICY 

NUMBER FAILED 



AVERAGE 

10.72 

10.70 

STANDARD DEVIATION 

2.82 

2.83 

REPAIR COST 



AVERAGE 

S 8.807. 

S 8,821. 

STANDARD DEVIATION 

S 2,3IS. 

S 2,298. 

PENALTY COST 



AVERAGE 

S 11,423. 

S 11.441. 

STANDARD DEVIATION 

S 5,591. 

S 5,582. 

TOTAL COST 



AVERAGE 

S 20.229. 

S 20.201. 

STAN DA R D D E V1 ATI O N 

S 7.146. 

S 7,116. 

PERCENT AVAILABIL¬ 
ITY 

99.9121 

99.9120 


Table 10. 20 YEAR STARTER SIMULATION RESULTS 
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WITH OUT POLICY 


WITH POLICY 


NUMBER FAILED 



AVERAGE 

4.79 

4.54 

STANDARD DEVIATION 

1.76 

1.90 

REPAIR COST 



AVERAGE 

S 2,567. 

S 2,891. 

STANDARD DEVIATION 

S 947. 

S 815. 

PENALT Y COST 



AVERAGE 

S 3,086. 

S 2,973. 

STANDARD DEVIATION 

S 4,463. 

S 4,395. 

TOTAL COST 



AVERAGE 

S 5,652. 

S 5,863. 

STANDARD DEVIATION 

S 4,779. 

S 4.681. 

PERCENT AVAILABIL¬ 
ITY 

99.9762 

99.7771 


Table 11. 20 YEAR LINK SIMULATION RESULTS 


D. TOWARDS SYSTEM LEVEL ANALYSIS AND SIMULATION 

Over the next two years data is going to be collected on the RC-IRON tanks and 
control group. Key components like the 48 candidates could be tabulated separately, 
aggregating the maintenance data for the remaining components. The aggregated fail¬ 
ures could be treated as a single component or as several components treated by 
catagory. These notional components could be analized using the techniques of this 
study. Because only components are renewed and not the system, simulation may be the 
way to explore alternatives. As an example of how this might work we shall take the 
six components who exibited Increasing failure rates and bring them together as an op¬ 
erating system. The program, Appendix I, is similar to the component simulation of the 
previous section. The difference is that we now have multiple components and a new 
process has been introduced to tie them together as a system. 
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A partial verification of this program was accomplished using queing theory. Five 
components were used in this test. Failures for each component were modeled as a 
special case of the Wiebull distribution, the Exponential, with mean failure time equal 
to five days. By chosing the Exponential distribution and keeping the repair times small, 
also Exponential with mean equal to 2.4 hours or 1/10 day, we expect the results to re¬ 
semble a M/M/1 queue. This is because we have a Poisson arrival process and because 
when all components are stopped when one fails we have a single server. The long run 
expected availability of the server is 90%. The simulated availability is 90.8% after 1,000 
repetitions. This shows that the simulation is working properly, 

The six component simulation was run 500 times. The results for mileage, number 
failures, and down time are located in Table 9. 



AVERAGE 

STANDARD DEVIATION 

MILEAGE 

4758 

10.49 

DOWN TIME (DAYS) 

1.51 

1.41 

NUMBER OF FAILURES 

4.07 

1.74 


Table 12. 2 YEAR SIMULATION RESULTS, SIX COMPONENT SYSTEM 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


An age replacement policy is not recommended for inclusion in the RC-IRON pro¬ 
gram based upon the 48 component NTC data set. As the program is currently estab¬ 
lished, the tanks in Germany would not have accumulated enough mileage for the three 
identified components (starter, grip, and link) to be replaced. There is a chance that the 
overall results may change if the Germany SDC data is analyzed. The environment of 
operation is different in Germany than at NTC so the components may exhibit different 
life distributions. If the Army leadership was to assign more value to availability than 
even the stand-by penalty, components may merit replacement. If the components were 
to be replaced in the field under an age replacement policy for the entire life cycle of the 
tank, the link would be the only contender. It should be noted that the assumptions of 
this study should be reexamined even for this component. This should be done from a 
engineering stand point. The link is used in conjuction with the tank track, so other 
factors may explain its failure distribution. Even with no parts being recommended for 
this type of replacement policy, this study should be of value. The life distributions have 
been examined and the components may be ranked in several ways. This may provide 
input into the inspection process, in that components which merit increased attention 
have been identified. A nonparainctric analysis along the lines of this thesis would be 
beneficial if the orginal NTC data set of time on test was expanded. 

A more general problem is the problem of determining which tanks should be sub¬ 
jected to RC-IRON. By simulating the tanks operation with components of interest, a 
near optimum solution for the mileage of RC-IRON application may be determined 
This simulation could use the techniques and information of this study. It will also re- 








quire the analysis of SDC Germany data and future NTC Hardware Test results. This 
system simulation may be built upon the simplistic simulation of the previous chapter. 
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APPENDIX A. Ml SYSTEM DESCRIPTION 


A. MISSION OF THE ABRAMS TANK 

The mission of the Abrams tank system is to close with and destroy enemy forces 
by use of firepower, maneuver and shock effect. The Abrams tank, organic to armored 
battalions and armored cavalry squadrons, will normally operate as part of combined 
arms team of armor, infantry and artillery to accomplish this mission. [Ref. 9: p.2] 

B. MODELS 

1. MI 

This is the basic model of the tank. Chrysler Corporation was awarded Full 
Scale Engineering Development in 1976 and sold its tank building subsidiary’ to General 
Dynamics in 1982. By 1985, the end of production, 2,374 were made. 

2. Improved M1 (IPM1) 

This is an Ml with improved armour protection. A tc.,al of 894 were built from 
1984 to 19S6. 

3. M1A1 

A number of improvements were made to the IPM1 for this tank. These in¬ 
clude. gun (see firepower next section), crew environment control, suspension, and 
transmission. Deliveries of this tank began in 1987 and are scheduled through 1991. 
By the end of 1989, 2,330 had been produced. [Ref. 10: pp.742-751] 

C CAPABILITIES 
1. Firepower 

The Ml version of the Abrams tank is armed with a 105 millimeter rifled can¬ 
non, the L6S. This cannon is combat-proven and arms the tanks of many allies. The 







M1A1 version, which entered production in 1985, is armed with the M256 120 millimeter 
smoothbore cannon, an improvement of the German 120mm cannon. Equipped with 
the M829 armor-piercing, fin-stabilized, discarding sabot round, this cannon can pene¬ 
trate any known main-battle tank armor currently fielded. A multipurpose high- 
explosive anti-tank round is also carried. A digital fire control computer, coupled with 
a laser range finder, thermal sights and a turret stabilization system enable the Abrams 
to engage targets under all weather conditions and on the move. The tank is also armed 
with a .50 caliber commander's maclnnegun and two 7.62 mm machineguns. 

2 . Mobility 

The Abrams was the worlds' first fielded tank to be equipped with a gas turbine 
engine. This engine develops 1500 horsepower and is coupled to a hydraulic trans¬ 
mission. An advanced suspension system featuring rotary shock absorbers enables the 
Abrams to operate at a maximum governed cross-country' speed of 42 miles per hour. 

3. Survivability 

The highest priority in the design of the Abrams was the protection of the 
crew. Compartmentation of fuel and ammunition, nuclear,'chcmical/biological pro¬ 
tection and halon fire suppression systems have been incorporated. Improved armor, 
responsive speed and agility, grenade and engine smoke generators and a low silhouette 
all contribute to the survivability of the system. 

4 . Communications 

Crew intercommunications are provided by the AN/VIC-1 intercom system. 
Tactical radio communications are provided by the AN/VRC-12 family of radios, with 
a maximum two net capability. Provisions are being made for the additions of position 
navigation and digital communication in future models. 
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5. Maintenance 


Maintenance considerations played a key part in the design of the Abrams. 
Ease of power pack removal and installation is the primary example. Additionally, most 
other major components are designed for easy removal and installation after fault iso¬ 
lation by built-in test equipment (BITE) or by the standard test equipment-Ml 
(STE-M1). [Ref. 9: p.3] 
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APPENDIX B. MAINTENANCE LEVELS 


A. UNIT MAINTENANCE 

This is the lowest level, it includes maintenance task performed by operator,, crew, 
and unit personnel. It may be equated to the maintenance performed by the owner of 
a car and service station. Preventive checks and services to detect potential problem is 
a key component. Replacements are limited to small components which are quickly and 
easily replaced. 

B. INTERMEDIATE MAINTENANCE 

Entire units are devoted to this level of maintenance. The intermediate is further 
divided into direct support and general support. Direct support is provided on a repair 
and return basis to units that experience failures beyond their capability to repair. 
General support units rebuild components in support of the Army supply system. 

C. DEPOT MAINTENANCE 

This is the highest level and is performed at large fixed depot facilities. They provide 
rebuild and overhaul for both systems and components. 
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APPENDIX C. RELIABILITY FITS 

A. STARTER 


STARTER 

WEIBULL CUMULATIVE DISTRIBUTION FUNCTION, N=113 



0 5000 10000 15000 20000 

MILES 


Figure 12. Starter Distribution Fit 









STARTER 

WEIBULL PR0E3ABIUTY PLOT, N=113 



MILES 


Figure 13. Starter Percentile Plot 








ANALYSIS OF WEIBULL DISTRIBUTION FIT 

DATA : STARTER 

SELECTION : ALL 
X AXIS LABEL: MILES 
SAMPLE SIZE : 113 

CENSORING : GROUPED DATA (CENSORING IS IMPLICIT) 
FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 

CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 


CONF. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER 

ESTIMATE 

LOWER 

UPPER 

C 

0 

C (SHAPE) 

1. 1999 

0.91873 

1.4811 

0. 020572 

2. 9477E1 

0 (SCALE) 

4689. 7 

3580. 5 

5798. 8 

29.477 

3. 2013E5 




SAMPLE* 

FITTED 

GOODNESS OF FIT 

MEAN 

: 

1974. 4 

4411.4 

NOT AVAILABLE 

STD DEV : 

1280. 8 

3692.2 


SKEWNESS: 

0. 54258 

1.5213 


KURTOSIS: 

2.6874 

6.2366 


* BASED ON 

MIDPOINTS OF 

FINITE INTERVALS 


PERCENTILES 

SAMPLE* 

FITTED 


5 


471 

394. 56 


10 


1C .'.5 

718. 86 


25 


1660 

1660.4 


50 


3352.5 

3455.3 


75 


5301 

6156.9 


90 



9397. 5 


95 



11702 



* BASED ON TURNBULL'S ESTIMATE 
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l. TRANSMISSION 



Figure 14. Transmission Distribution Fit 











PERCENTILE 


TRANSMISSION 

WEIBULl PROBABILITY PLOT. N=91 



Figure 15. Transmissionr Percentile Plot 
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ANALYSIS OF WEIBULL DISTRIBUTION FIT 

DATA : TRANSMISSION 

SELECTION : ALL 
X AXIS LABEL: MILES 
SAMPLE SIZE : 91 

CENSORING : GROUPED DATA (CENSORING IS IMPLICIT) 

FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 

CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 

CONF. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 

PARAMETER ESTIMATE LOWER UPPER C 0 

C (SHAPE) 1.4455 0.9569 1.934 0.062105 2. 0929E2 

0 (SCALE) 7479.5 4952.5 10007 209.29 1. 6616E6 


SAMPLE* FITTED GOODNESS OF FIT 

MEAN : 2441.2 6784.8 NOT AVAILABLE 

STD DEV : 1325.1 4767 

SKEWNESS: 0.040961 1.1367 

KURTCSIS: . 2.2564 4.6203 

* BASED ON MIDPOINTS OF FINITE INTERVALS 

PERCENTILES SAMPLE* FITTED 

5: 665 956.21 

10: 1777.5 1576.7 

25: 3153.5 3158.9 

50: 5122 5804.3 

75: 9375.9 

90: 13319 

95: 15978 

* BASED ON TURNBULL’S ESTIMATE 
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C. GRIP 


GRIP 

WEIBUIL CUMULATIVE DISTRIBUTION FUNCTION, N«118 



Figure 16. Grip Distribution Fit 
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PERCENTILE 


GRIP 

WEIBULL PROBABILITY PLOT, N=118 



TOO 1000 10000 

MILES 


Figure 17. Grip Percentile Plot 
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ANALYSIS OF WEIBULL DISTRIBUTION FI1 


DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORING 
FREQUENCIES 
EST. METHOD 
CONF METHOD 


GRIP 

ALL 

MILES 

118 

GROUPED DATA (CENSORING IS IMPLICIT) 
1 

MAXIMUM LIKELIHOOD 
ASYMPTOTIC NORMAL APPROXIMATION 




CONF. INTERVALS 

COVARIANCE 

MATRIX OF 



(95 PERCENT) 

PARAMETER ] 

ESTIMATES 

PARAMETER 

ESTIMATE 

LOWER UPPER 

C 

0 

C (SHAPE) 

1. 5539 

1.1948 1.913 

0.033552 

3.0428E1 

0 (SCALE) 

4369. 9 

3538.7 5201.1 

30. 428 

1.7976E5 



SAMPLE* 

FITTED 

GOODNESS OF FIT 

MEAN : 

2060.9 

3929.2 

NOT AVAILABLE 

STD DEV : 

1072.3 

2582. 6 


SKEWNESS: 

0. 10384 

1. 0109 


KURTOSIS: 

2. 2689 

4.1929 


* BASED ON 

MIDPOINTS OF 

FINITE INTERVALS 


PERCENTILES 

SAMPLE* 

FITTED 


5: 

365. 5 

646.17 


10: 

1108.5 

1026.9 


25: 

2039 

1960 


50: 

3160 

3451. 7 


75: 


5392. 1 


90: 


7474.3 


95: 


8853. 6 


* BASED ON 

TURNBULL'S ESTIMATE 
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CUMULATIVE PROBABILITY 
O 2 D.4 0.6 


D. NOZZLE 


NOZZLE 

WEIBULL CUMULATIVE DISTRIBUTION FUNCTION. N=117 



Figure 18. Nozzle Distribution Fit 
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ANALYSIS OF WEIBULL DISTRIBUTION FIT 


DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORING 
FREQUENCIES 
EST. METHOD 
CONF METHOD 


NOZZLE 

ALL 

MILES 

117 

GROUPED DATA (CENSORING IS IMPLICIT) 
1 

MAXIMUM LIKELIHOOD 
ASYMPTOTIC NORMAL APPROXIMATION 


PARAMETER ESTIMATE 
C (SHAPE) 1.3384 

0 (SCALE) 4438.7 


CONF. INTERVALS 
(95 PERCENT) 

LOWER UPPER 

1.0287 1.6482 

3490.1 5387.2 


COVARIANCE MATRIX OF 
PARAMETER ESTIMATES 
C 0 

0.024967 2.9137E1 

29.137 2.3412E5 


SAMPLE* FITTED GOODNESS OF FIT 

MEAN : 1925.3 40/6.5 NOT AVAILABLE 

STD DEV : 1223.4 3077 

SKEWNESS: 0.91959 1.2864 

KURTOSIS: 3.4006 5.1839 

* BASED ON MIDPOINTS OF FINITE INTERVALS 

PERCENTILES SAMPLE* FITTED 
5: 571 482.47 

10: 1057 826.11 

25: 1489 1749. 7 

50: 3457 3375.4 

75: 5124 5665.5 

90: 5124 827/'. 2 

95: 5124 10076 

* BASED ON TURNBULL'S ESTIMATE 
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CUMULATIVE PROBABILITY 
0.2 0.4 D.6 0.8 


E. DISTBOX 


DISTBOX 

WEIBULL CUMULATIVE DISTRIBUTION FUNCTION, N=94 



Figure 20. Distbox Distribution Fit 
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DISTBOX 

WEiBULL PROBABILITY PLOT, N=94- 



100 1000 10000 

MILES 


Figure 21. Distbox Percentile Plot 
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ANALYSIS OF WEIBULL DISTRIBUTION FIT 


I DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORING 
FREQUENCIES 
EST. METHOD 
CONF METHOD 


DISTBOX 

ALL 

MILES 

94 

GROUPED DATA (CENSORING IS IMPLICIT) 
1 

MAXIMUM LIKELIHOOD 
ASYMPTOTIC NORMAL APPROXIMATION 




CONF. INTERVALS 

COVARIANCE MATRIX OF 



(95 PERCENT) 

PARAMETER ESTIMATES 

PARAMETER 

ESTIMATE 

LOWER UPPER 

C 0 

C (SHAPE) 

1. 5315 

1.0412 2.0217 

0.062543 1.5897E2 

0 (SCALE) 

6858. 6 

4828.2 8889.1 

158.97 1.0727E6 




SAMPLE* 

FITTED 

GOODNESS OF FIT 

MEAN 

: 

2705.2 

6176.7 

NOT AVAILABLE 

STD DEV : 

1535. 9 

4114.4 


SKEWNESS: 

0.0087015 

1. 0358 


KURTOSIS: 

1. 904 

4.2719 


* BASED ON 

MIDPOINTS OF 

FINITE INTERVALS 


PERCENTILES 

SAMPLE* 

FITTED 


5 


975. 5 

986. 16 


10 


1431. 5 

1577. 9 


25 


3633.5 

3040.4 


50 


5111.5 

5398. 9 


75 


5599 

8489. 2 


90 



11824 


95 



14040 


* BASED ON 

TURNBULL'S ESTIMATE 



£ 
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CUMULATIVE PROBABILITT 
D.2 0.4 0.6 0.3 


F. LINK 


LINK 

WEIBULL CUMULATIVE DISTRIBUTION FUNCTION, N=8B 



0 10000 20000 30000 

MILES 


Figure 22. Link Distribution Fit 
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PERCENTILE 



figure 23. Link Percentile Plot 
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ANALYSIS OF WEIBULL DISTRIBUTION FIT 


DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORINC 
FREQUENCIES 
EST. METHOD 
CONF METHOD 


LINK 

ALL 

MILES 

88 

GROUPED DATA (CENSORING IS IMPLICIT) 
1 

MAXIMUM LIKELIHOOD 
ASYMPTOTIC NORMAL APPROXIMATION 


PARAMETER ESTIMATE 
C (SHAPE) 1.2801 

0 (SCALE) 10239 


CONF. INTERVALS 
(95 PERCENT) 

LOWER UPPER 

0.78478 1.7754 

5049.7 15429 


COVARIANCE MATRIX OF 
PARAMETER ESTIMATES 
C # 

0.063837 5.0284E2 

502.84 7. 0077E6 


GOODNESS OF FI" 
NOT AVAILABLE 


PERCENTILES SAMPLE* FITTED 


5: 1196.5 1006 

10: 1753 1765.2 

25: 3592.5 3868.8 

50: 7689. 9 

75: 13215 

90: 196 *4 

95: 24127 


* BASED ON TURNBULL'S ESTIMATE 


SAMPLE* FITTED 
MEAN : 2071.2 9486.9 

STD DEV : 1137.5 7467 

SKEWNESS: 0.57276 1.3783 

KURTOSIS: 2.5741 5.5727 

* BASED ON MIDPOINTS OF FINITE INTERVALS 









APPENDIX D. DELAY FITS 


A. STARTER 


STARTER DELAY 

LOGNORMAL OENSITY FUNCTION. N=57 



Figure 24. Starter Delay Histogram 









STARTER DELAY 

LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=57 



l. 


Figure 25. Starter Delay Distribution Fit 
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PERCENTILE 


STARTER DELAY 

LOGNORMAL PROBABILITY PLOT, N=57 



0.1 0.2 0.4 12 4 

DAYS 


Tigure 26. Starter Delay Percentile Plot 
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ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA 

STARTER 



SELECTION 

ALL 



X AXIS LABEL 

DAYS 



SAMPLE SIZE 

57 



CENSORING 

NONE 



FREQUENCIES 

1 



EST. METHOD 

MAXIMUM 

LIKELIHOOD 


CONF METHOD 

EXACT 





CONF. INTERVALS 

COVARIANCE MATRIX OF 



(95 PERCENT) 

PARAMETER ESTIMATES 

PARAMETER ESTIMATE LOWER UPPER 

MU SIGMA 

MU 

0.59102 

0.80521 0.37683 

0.011212 0 

SIGMA 

0.79943 

0.6809 0.98954 

0 0.0056061 



SAMPLE 

FITTED 

GOODNESS OF FIT 

MEAN : 

0.79685 

0.76225 



STD DEV : 

0.94387 

0.72103 

KOLM-SMIRN 

0.18602 

SKEWNESS: 

3. 6941 

3.6841 

SIGNIF 

0.038716 

KURTOSIS: 

18.884 

34.264 

CRAMER-V M 

0.34074 




SIGNIF 

< . 15 

PERCENTILES 

SAMPLE 

FITTED 

ANDER-DARI. 

1. 5772 

5: 

0. 125 

0.14864 

SIGNIF 

> . 15 

10 : 

0 . 2 

0.19876 



25: 

0. 4 

0.32304 

KS, AD, AND CV SIGNIF. LEVELS NOT 

50: 

0. 5 

0.55376 

EXACT WITH ESTIMATED PARAMETERS. 

75: 

0. 75 

0.94928 



90: 

2 

1.5429 



95: 

2. 25 

2.0631 
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RELATIVE FREQUENCY 
0.2 0.4 


B. TRANSMISSION 


TRANSMISSION DELAY 
LOGNORMAL DENSITY FUNCTION, N«29 



Figure 27. Transmission Delay Histogram 


64 










‘UMULATTVE PROBABILITY 



TRANSMISSION DELAY 

LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=29 



Figure 28. Transmission Delay Distribution fit 
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PERCENTILE 


TRANSMISSION DELAY 
LOGNORMAL PROBABILITY PLOT, N=29 



Figure 29. Transmission Delay Percentile Plot 
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ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORING 
FREQUENCIES 
EST. METHOD 
CONF METHOD 


TRANSMISSION 

ALL 

DAYS 

29 

NONE 

1 

MAXIMUM LIKELIHOOD 
EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER 

ESTIMATE 

LOWER 

UPPER 

MU 

SIGMA 

MU 

0.15064 

0. 51526 

0. 21398 

0.030588 

0 

SIGMA 

0. 94134 

0. 76064 

1. 2964 

0 

0.015294 



SAMPLE 

FITTED 

GOODNESS OF FIT 

MEAN : 

1. 295 

1.3403 



STD DEV : 

1. 2194 

1. 6016 

CHI-SQUARE 

0.23296 

SKEWNESS: 

1.5456 

5.2913 

DEG FREED 

1 

KURTOSIS: 

4.7651 

78.062 

SIGNIF 

0.62934 




KOLM-SMIRN 

0.13149 

PERCENTILES 

SAMPLE 

FITTED 

SIGNIF 

0. 69773 

5: 

0. 2 

0.18265 

CRAMER-V M 

0. 054591 

10: 

0. 25 

0.25722 

SIGNIF 

> . 15 

25: 

0. 5 

0.45584 

ANDF.R-DARL 

0.29512 

50: 

1 

0.86016 

SIGNIF 

> . 15 

75: 

1. 5652 

1. 6231 



90: 

3. 5 

2.8764 

KS, AD, AND CV SIGNIF. LEVELS NOT 

95: 

4 

4.0506 

EXACT WITH ESTIMATED PARAMETERS. 



CHI-SQUARE 

GOODNESS OF 

FIT TABLE 


LOWER 

UPPER 

OBS 

EXP 

0-E 

((0-E)*2)*E 

-INF. 

0. 61746 

13 

13.875 

0.87491 

0.05517 

0. 81746 

1.6349 

9 

7.9432 

1. 0568 

0.1406 

1.6349 

2.4524 

3 

3.3254 

0.32541 

0. 031842 

2.4524 

+ INF. 

4 

3. 8565 

0.14351 

0. 0053407 

TOTAL 


29 

29 


0.23296 
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RELATIVE FREQUENCY 


C. GRIP 


GRIP DELAY 

LOGNORMAL DENSITY FUNCTION, N=5< 



Figure 30. Grip Delay Histogram 
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CUMULATIVE PROBABILITY 


GRIP DELA' 



Figure 31. Grip Delay Distribution Fit 
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Figure 32. Grip Delay Percentile Plot 
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ANALYSIS Cl LOGNORMAL DISTRIBUTION FIT 

DA"'A : GRIP 

SELECT'ON : ALL 
X AXIS LABEL: DAYS 
SAMPLE SIZE : 54 

CENSORING : NONE 
FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 

COMF. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 

PARAMETER ESTIMATE LOWER UPPER MU SIGMA 


MU 

0.60291 

0.84539 

0.36042 

0. 014336 

0 

SIGMA 

0.87986 

0.74655 

1. 0965 

0 

0. 007168 




SAMPLE 

FITTED 

GOODNESS OF 

FIT 

MEAN 

: 

1. 0154 

0.80587 



STD DEV : 

2. 5633 

0. 87121 

KOLM-SMIRN 

0.17048 

SKEWNESS: 

6. 588 

4. 5068 

SIGNIF 

0. 086673 

KURTOSIS: 

46. 611 

53.634 

CRAMER-V M 

0. 25423 





SIGNIF 

> . 15 

PERCENTILES 

SAMPLE 

FITTED 

ANDER-DARL 

1. 3338 

5 


0. 125 

0.12868 

SIGNIF 

> . 15 

10 


0.15385 

0.1771/ 



25 


0.33333 

0.30237 

KS, AD, AND 

IV SIGNIF. LEVELS NOT 

50 


0. 5 

0.54722 

EXACT WITH ESTIMATED PARAMETERS. 

75 


1 

0. 99033 



90 


1. 2 

1. 6902 



95 


2 

2. 3272 
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RELATIVE FREQUENCY 


D. NOZZLE 


NOZZLE DEWY 

LOGNORMAL DENSiTFUNCTION, N=5B 



Figure 33. Nozzle Delay Histogram 









CUMULATIVE PROBABILFTf 
0.2 0.4 0.6 


NOZZLE DELAY 

LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N»56 



I 


Figure 34. Nozzle Del*> Distribution Fit 








PERCENTILE 


NOZZLE DELAY 

LOGNORMAL PROBABILITY PLOT, N=56 



DAYS 


Figure 35. Nozzle Delay Percentile Plot 


7J 


























ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA : NOZZLE 

SELECTION : ALL 
X AXIS LABEL: DAYS 
SAMPLE SIZE : 56 

CENSORING : NONE 
FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 




CONF. INTERVALS 

COVARIANCE MATRIX OF 



(95 PERCENT) 

PARAMETER ESTIMATES 

PARAMETER 

ESTIMATE 

LOWER UPPER 

MU SIGMA 

MU 

0. 69437 

0.876 0.51275 

0.0080534 0 

SIGMA 

0.67156 

0.57128 0.83309 

0 0.0040267 


MEAN 

. 

SAMPLE 

0.61794 

FITTED 

0. 6257 

GOODNESS OF FIT 

STD DEV : 

0. 42191 

0.47234 

CHI-SQUARE : 11.675 

SKEWNESS: 

1, 4959 

2.6949 

DEG FREED: 3 

KURTOSIS: 

5.2464 

18.205 

SIGNIF : 0.0085837 

KOLM-SMIRN : 0.19716 

PERCENTILES 

SAMPLE 

FITTED 

SIGNIF : 0.025722 

5 


0. 13208 

0.16543 

CRAMER-V M : 0.34086 

10 


0.23077 

0. 21116 

SIGNIF : < .15 

25 


0.31667 

0.31755 

ANDER-DARL : 1.552 

50 

75 


0. 5 

0.92857 

0.49939 

0.78536 

SIGNIF : > .15 

90 


1.1429 

1. 181 

KS, AD, AND CV SIGNIF. LEVELS NOT 

95 


1. 5 

1. 5075 

EXACT WITH ESTIMATED PARAMETERS. 



CHI-SQUARE 

GOODNESS OF 

FIT TABLE 


LOWER 

UPPER 

OBS 

EXP 

0-E 

((0-E)*2)*E 

-INF. 

0.27381 

12 

10.384 

1. 6158 

0.25143 

0.27381 

0.54762 

27 

20. 673 

6.3265 

1. 9361 

0.54762 

0.82143 

2 

12. 1 

10. 1 

8.4306 

0 8214j 

1.0952 

8 

6. 06 

1. 94 

0.62102 

1.0952 

1. 3691 

4 

3. 0536 

0. 94643 

0. 29334 

1. 3691 

+INF. 

3 

3. 7287 

0.72874 

0. 14242 

TOTAL 


56 

56 


11.675 









RELATIVE FREQUENCY 


E. Dio7?OX 


DISTBOX DELAY 

LOGNORMAL DENSITY FUNCTION, N=32 



0 2 4 6 

DATS 


Figu-e 36. Distbox Deiay Histogram 
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Figure 37. Distbox Delay Distribution Fit 
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Figure 38. Distbox Delay Percentile Plot 


i 
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ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA : DISTBOX 

SELECTION : ALL 
X AXIS LABEL: DAYS 
SAMPLE SIZE : 32 

OEN&UKING : NONE 
FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER 

ESTIMATE 

LOWER 

UPPER 

MU 

SIGMA 

MU 

0.57369 

0. 916 

0.23139 

0.027274 

0 

SIGMA 

0.93421 

0.76092 

1. 2624 

0 

0.013637 




SAMPLE 

FITTED 

GOODNESS OF 

FIT 

MEAN 

: 

0.90017 

0. 8717 



STD DEV : 

1.1137 

1. 029 

KOLM-SMIRN 

0.14462 

SKEWNESS: 

2.9283 

5.1864 

SIGNIF 

0.51498 

KURTOS 

IS: 

12.181 

74.429 

CRAMER-V M 

0.092307 





SIGNIF 

> . 15 

PERCENTILES 

SAMPLE 

FITTED 

ANDER-DARL 

0. 47943 

5 


0. 13766 

0.12116 

SIGNIF 

> . 15 

10 


0. 154 

0.17015 



25 


0. 275 

0.30013 

KS, AD, AND 

7V SIGNIF. LEVELS NOT 

50 


0. 5 

0. 56344 

EXACT WITH ESTIMATED PARAMETERS. 

75 


1 

1.0577 



90 


2 

1.8658 



95 


3.4286 

2. 6203 
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RELATIVE FREQUENCY 


F. LINK 


LINK DELAY 

LOGNORMAL DENSITY FUNCTION, N=22 



Figure 39. Link Delay Histogram 
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CUMULATIVE PROBABILITY 
0.2 D.4 0.6 0.8 


LINK DELAY 


LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=22 



J 


Figure 40. Link Delay Distribution Fit 


8! 







LINK DELAY 

LOGNORMAL PROBABILITY PLOT, N=22 


Figure 4i. Link Delay Percentile Plot 
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ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORING 
FREQUENCIES 
EST. METHOD 
CONF METHOD 


LINK 

ALL 

DAYS 

22 

NONE 

1 

MAXIMUM LIKELIHOOD 
EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER 

ESTIMATE 

LOWER 

UPFER 

MU 

SIGMA 

MU 

0. 351 

0.88451 

0.1825 

0.062811 

0 

SIGMA 

1.1755 

0. 9256'' 

1. 7195 

0 

0.031405 




SAMPLE 

FITTED 

GOODNESS OF FIT 

MEAN 

: 

2. 091 

1.h048 



STD DLV : 

5.6498 

2. 426 

KOLM-SMIRN 

0.25086 

SKEWNESS: 

4.1585 

10.331 

SIGNIF 

0. 12542 

KURTOSIS: 

18.856 

422. 36 

CRAMER-V M 

0. 23969 





SIGNIF 

> . 15 

PERCENTILES 

SAMPLE 

FITTED 

ANDER-DARL 

1. 2127 

5 


0. 25 

0.10178 

SIGNIF 

> . 15 

10 


0. 25 

C. 15o83 



25 


0. 4 

0.316o9 

KS, AD, AND CV SIGNIF. LEVELS NOT 

50 


0. 5 

0. 70398 

EXACT WITH ESTIMATED PARAMETERS. 

75 


1 

1.5551 



90 


3 

3. 1762 



95 


4 

4. ft/;- 
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RELATIVE FREQUENCY 


APPENDIX E. REPAIR TIME UTS 



Starter keuair Histogram 







CUMULATIVE PROBABILfTT 


STARTER REPAIR TIME 

LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=57 



HOURS 


Figure 43. Starter Repair Distribution Fit 
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Figure 44. Starter Repair Percentile Plot 



















ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA STARTER 

SELECTION : ALL 
X AXIS LABEL: HOURS 
SAMPLE SIZE : 57 

CENSORING : NONE 
FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER 

ESTIMATE 

LOWER 

UPPER 

MU 

SIGMA 

MU 

0. 55501 

0. 39739 

0. 71262 

0.0060716 


SIGMA 

0.58829 

0. 50106 

0. 72818 

0 

0. 0030358 


SAMPLE 

MEAN : 2.0018 

FITTED 

2. 071 

GOODNESS OF 

FIT 

STD DEV : 1.0302 

1. 3318 

CHI-SQUARE 

7. 5867 

SKEWNESS: 1.3499 

2. 1951 

DEG FREED 

3 

KTRTOSIS: 5.7311 

12.635 

SIGNIF 

KOLM-SMIRN 

0.055365 

0. 1791 

PERCENTILES SAMPLE 

FITTED 

SIGNIF 

0. 051622 

5: 0. 9 

0.66176 

CRAMER-V M 

0.33067 

10 : 1 

0.81954 

SIGNIF 

< . 15 

25: 1. 3 

1. 1716 

ANDER-DARL 

1. 9312 

50: 2 

75: 2 

1. 742 

2. 5899 

SIGNIF 

< . 15 

90: 3 

3. 702o 

KS, AD, AND 

IV SIGNIF. LEVELS NOT 

95: 4 

4.565 3 

EXACT WITH ESTIMATED PARAMETERS. 


CHI-SQUARE GOODNESS OF FIT TABLE 


LOVER 

UPPER 

OBS 

EXP 

0-E 


-INF. 

0 84286 

9 

*. 

6.1899 

4. 1899 

2. 63e 1 

0.84286 

1. t> 5 5 7 

: i 

21 . 0-2 

0. u-*2 19 7 

C. 050Oft,cIS 

1. 1 6 5 7 

2 . 52St? 

-? y 

4. A 

1 , t?-- 

6.2357 

2 .633' 

2.5266 

3. 3 7!-* 

8 

7. 5 4 f> 1 

0. 45 387 

0 -'299 

3. 3714 

4 2L3 

** 

3. 6*2 5 

0 33747 

C. C.J l' 4 *-. 

4. 2143 

TOTAL 

+ INF. 

c 7 

3. 7v-9 

5' 

2 . 79~»^ 

2. (. 5 5- 

7 55c’ 


S' 








RELATIVE FREQUENCY 

0.2 0.3 D.4 0.5 


B. LINK 


LINK REPAIR TIME 

LOGNORMAL DENSITY FUNCTION, N=22 



HOURS 


Figure 45. Link Repair Histogram 
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CUMULATIVE PROBABILITY 
O.Z 0.4 0.6 


LINK REPAIR TIME 

LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=22 



Figure 46. Link Repair Distribution Fit 
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PERCENTILE 


LINK REPAIR TIME 
LOGNORMAL PROBABILITY PLOT, N=22 



Figure 47. Link Repair Percentile Plot 
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ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 






DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORING 
FREQUENCIES 
EST. METHOD 
CONF METHOD 


LINK 

ALL 

HOURS 

22 

NONE 

1 

MAXIMUM LIKELIHOOD 
EXACT 


CONI'. INTERVALS COVARIANCE MATRIX OF 

(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER 

ESTIMATE 

LOWER 

UPPER 

MU 

SIGMA 

MU 

0.96595 

0. 67686 

1. 253 

0.018188 

0 

SIGMA 

0.63257 

0. 49812 

0. 92531 

0 

0.0090941 




SAMPLE 

FITTED 

GOODNESS OF FIT 

MEAN 

: 

3.2455 

3. 2092 



STD DEV : 

2.4705 

2. 2511 

CHI-SQUARE 

1. 2103 

SKEWNESS: 

1. 9599 

2.4495 

DEG FREED 

1 

KURTOSIS: 

6.4614 

15.277 

SIGNJ.F 

0.27117 





KOLM-SMIRN 

0.16686 

PERCENTILES 

SAMPLE 

FITTED 

SIGNIF 

0.57263 

5 


1 

0.92796 

CRAMER-V M 

0.079711 

10 


1.5 

1.1679 

SIGNIF 

> . 15 

25 


2 

1. 7151 

ANDER-DARL 

0.47637 

50 


2. 25 

2.6273 

SIGNIF 

> . 15 

75 


4 

4. 0246 
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4. 5 

5.9105 

KS, AD, AND CV SIONIF. LEVELS NOT 

95 


9 

7 4384 

EXACT WITH ESTIMATED PARAMETERS. 



CHI-SQUARE 

GOODNESS 

OF FIT TABLE 

LOWER 

UPPER 

OBS 

EXP 

0-E 

((0-E)*2)*E 

-INF. 

1.7167 

5 

5.5122 

0.51216 

0.047566 

1. 7167 

3.4333 

9 

9.0927 

0.092724 

0.00094556 

3. 4333 

5. 15 

6 

4.2345 

1.7655 

0.7361 

5. 15 

+INF. 

2 

3.1606 

1. 1606 

0 4262 

TOTAL 


22 

22 


1.2108 
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APPENDIX F. PROGRAM USED TO DETERMINE AGE 


REPLACEMENT AND COST 


Code: APL 

Programmer: O. Uyar, provided by author of [Ref. 11] 
Date: August 1990 


V SIM;CliC2iIiJiFXiXiTiXAiYA',CiXMIN;YMIN 
-1- fl THIS PROGRAM SIMULATES THE COST FUNCTION (EQUATION 2.4) TO FIND 
-2- A MINIMUM VALUE (YMIN ) OF THE COST FUNCTION AND CORRESPONDING AGE , 
-3- a REPLACEMENT TIME ( XMIN ) FOR THAT POINT . AFTER FINDING MINIMUM 

[4] r PAZRFS INSIDE THE LOOP1 IT REPEAT THE PROCEDURE 30C TIMES INSIDE 

[5] A THE L00P2. FINALLY, THE PROGRAM GIVES THE AVERAGE VALUES FOR 

[6] r BOTH MINIMUM POINT AS AXST AND ACST . 

[7] T+ (t500Q)*100 

-8- fl THIS GIVES US A VECTOR OF T(0.01, 0.02, ..., 50) TO CALCULATE 
[93 R FIRST £7(0.01) AND THEN £7(0.02) UP TO £7(50) OF 5000 COST VECTOR. 

[10] fl INITIALIZATION ... 

[11] fl UNPLANNED AND PLANNED REPLACEMENT COST MUST EE GIVEN BY THE USER. 

[12] £71-«-5 

[13] £72-*-1 
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[14] XA+iO 

[15] YA +lO 

[16] J +0 

[17] fi J IS THE INCREMENT OF THE LOOP2 J= 1, 2, ...,300 

[18] o MODEL ... 

[19] L00P21 

[20] X+-5000 WEIRAND 2 2.2567587 

[21] p LINE 14, GENERATES 5000 SYSTEM LIFETIMES FROM 

[22] P WEI ( ALPHA- 2.0 , BETA-2 .2567587 )AS VECTOR X. HERE BETA VALUE 

[23] p REPRESENTS 1 OVER LAMBDA =( 1+0.44311346 ). 

[24] p FOi? GAMMA DISTRIBUTION LINE 14 CAW BE SWITCH WITH 

[25] X*5C00 WEIRAND 4 0.5 FOR GAMMA ( P=4 , THETA-O.b ). 

[26] J+J+l 

[27] C+iO 

[28] I«-0 

[29] p I IS THE INCREMENT OF THE INNER LOOP 1=1 , 2, ...,5000 

[30] LOOP !: 

[21] I«-J+l 

[32] R C 15 T5F SIMULATED COST FUNCTION 

[33] C+C, (((C2*(i-FX)) + (Clx(FX«-(( + /X£r[I] )* 5000 ))))+((+/(XLW] )) 
-34- p IN THE FIRST LOOP C VECTORS OBTAIN FOR EACH T 

[35] »(I<5000 VLOOP1 

[36] IMIN+l/C 

[37] XMIN-^T 11 *kCl 
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[38] A YMIN i THE MINIMUM VALUE OF THE COST FUNCTION FOR SPECIFIC T 

[39] A XMIN : THE CORRESPONDING AGE REPLACEMENT TIME (T) 

[40] XA+XA,XMIN 

[41] Y A-*-Y A, YMIN 

[42] a XA : THE VECTOR OF THE ACE REPLACEMENT TIMES (300) 

[43] A YA : THE VECTOR OF THE YMIN (300) 

[44] -*-(«7<300 )/L00P2 

( 

[45] AXST+(+/XA)*pXA 

[46] ACST+(+/YA)*pYA 

-47- a AXST : THE AVERAGE VALUE OF THE AGE REPLACEMENT TIMES AFTER 300 REP 

-48- a ACST i THE AVERAGE VALUE OF THE YMIN AFTER 300 REPEATITIONS . 

V 
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APPENDIX G. SCALED TIME ON TEST PLOT PROCEDURE AND APL 


CODE 


Code: APL 

Programmer: J. Wilhelm 
Date: August 1990 


STEPS 

1) MAKE A 101 UNIT VECTOR FROM 0 TO MAXIMUM 
MILAGE WITH INTERVALS H. 

» XX-*- 0,ffx(xi00) 

2) USE THE XX VECTOR TO FIND A WEIBULL CDF 

* 

VECTOR. 

P+SHAPE SCALE WEICDF XX 

3) CALCULATE THE SURVIVAL VECTOR 
F3AR+ 1 - P 

STARTER INPUTS 

H<r 200 

XX* 0,200*(il00) 

P+1.1999 4689.7 WEICDF XX 
FEAR+ 1 - P 
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LINK INPUTS 


H* 300 

XX* 0,300x(xl00) 

P+1.2801 10239. WEICDF XX 
FBAR* 1 - P 


4) RUN THIS PROGRAM TO INTEGRATE THE SURVIVOR 
FUNCTION FROM 0 TO EACH POINT IN THE VECTOR. 
THE RESULT IS A VECTOR K. 

VBLDG [□] V 

V BLDG 

[1] J+2 

C2] K*0 

[3] LOOP:+ENDLOOP IF J>101 

[4] PART*I+FBAR 

[5] J*H SIMPSON PART 

[6] K*K,J 

[7] 1*1+1 

[8] *LOOP 

[9] ENDLOOP: 

[10] 

V 
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5 ) STANDARDIZE K WITH THE MEAN. 

STD-*- K * 4411.4 

RESULTS 

STARTER 

STD-*- K * 4411.4 
LINK 

STD-*r K * 9486.9 

6 ) PLOT STD VERSUS P TO OBTAIN THE 
SCALED TIME ON TEST. 
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APPENDIX H. COMPONENT SIMULATION 


Code: SIMSCR1PT II.5 
Programmer: J, Wilhelm 
Date: August 1990 


’' PART RAM MODEL 

PREAMBLE 

NORMALLY, MODE IS UNDEFINED 
DEFINE .MILES TO MEAN MINUTES 
PROCESSES INCLUDE STOP. SIM AND PART 
DEFINE FAIL.SHAPE AS A REAL VARIABLE 
DEFINE FAIL. SCALE AS A REAL VARIABLE 
DEFINE DELAY. MU AS A REAL VARIABLE 
DEFINE DELAY.SIG AS A REAL VARIABLE 
DEFINE REP. MU AS A REAL VARIABLE 
DEFINE REP.SIG AS A REAL VARIABLE 
DEFINE SEED1 AND SEED2 AS INTEGER VARIABLES 
DEFINE SEED3 AS A INTEGER VARIABLE 
DEFINE NO. OF.PARTS AND NO. FAIL AS INTEGER VARIABLES 
DEFINE REPL. NO AS A INTEGER VARIABLE 
DEFINE N AS A INTEGER VARIABLE 
DEFINE NOW. WORK AS A INTEGER VARIABLE 
DEFINE PRCNT. UP AS A REAL VARIABLE 
DEFINE T. FAIL AS A INTEGER VARIABLE 
DEFINE RUN. LENGTH AS A REAL VARIABLE 
DEFINE MILES.T AS A REAL VARIABLE 
DEFINE MILES.N AS A REAL VARIABLE 
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DEFINE DOWN. T AS A REAL VARIABLE 
DEFINE DOWN. N AS A REAL VARIABLE 
DEFINE PART. COST AS A REAL VARIABLE 
DEFINE STAND. BY. COST AS A REAL VARIABLE 
DEFINE MN. HR. COST AS A REAL VARIABLE 
DEFINE PEN.COST AS A REAL VARIABLE 
DEFINE REP. COST AS A REAL VARIABLE- 
DEFINE TOT. COST AS A REAL VARIABLE 
DEFINE D. COST AS A REAL VARIABLE 
DEFINE R. COST AS A REAL VARIABLE 
DEFINE T. COST AS A REAL VARIABLE 
TALLY AVG. MILES. T AS THE MEAN, 

SIG.MILES. T AS THE STD. DEV OF MILES. T 
TALLY AVG. DOWN.T AS THE MEAN, 

SIG.DOWN.T AS THE STD.DLV OF DOWN.T 
TALLY AVG.T. FAIL AS THE MEAN, 

SIG. T. FAIL AS THE STD. DEV OF T. FAIL 
TALLY AVG. R. COST AS THE MEAN, 

SIG.R. COST AS THE STD. DEV OF R.COST 
TALLY AVG. D. COST AS THE MEAN, 

SIG.D.COST AS THE STD. DEV OF D.COST 
TALLY AVG. T. COST AS THE MEAN, 

SIG. T. COST AS THE STD. DEV OF T. COST 
ACCUMULATE AVAIL AS THE AVERAGE OF NOW. WORK 


MAIN 

LET MINUTES. V = 6.5373 / 24. ''MILES 
CALL READ. DATA 
FOR N = 1 TO 5000,DO 
CALL INITIALIZE 

ACTIVATE A STOP. STM IN RUN. LENGTH DAYS 




START SIMULATION 




LOOP 


END 

ROUTINE INITIALIZE 
" STARTER 
*' CREATE A PART 
'' LET FAIL. SHAPE - 1.1999 
" LET FAIL. SCALE = 4689.7 
" LET DELAY. MU = .59102 
'' LET DELAY. SIG = . 79943 
" LET REP.MU =.55501 
'' LET REP.SIG = .58829 

'' LET SEED1 = 8 

'' LET SEED2 = 2 

'' LET SEED3 = 4 

'' PART.COST = 794. 

1 1 ********* 

'' LINK 
CREATE A PART 
LET FAIL. SHAPE = 1.2801 
LET FAIL. SCALE = 10239.0 
LET DELAY. MU = .3510 
LET DELAY.SIG = 1. 1755 
LET REP. MU =.96595 
LET RF.P. SIG = . 6325 7 

LET SEEDi = 9 

LET SEED2 = 1 

LET SEED3 = 5 

PART. COST = 488. 

''****** 

ACTIVATE THIS PART NOW 


END 
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ROUTINE READ. DATA 

NO.OF. PARTS = 1 
'' PRINT 1 LINE THUS 
" HOW LONG SHALL WE RUN (IN DAYS)? 

'' READ RUN. LENGTH 

RUN. LENGTH = 730. 0 " 730 DAYS IN TWO YEAR OR 7300 IN 20 

STAND. BY. COST = 1781. "PENALTY COST PER DAY 
MN. HR. COST = 50. " COST PER REPAIR HOUR 

END 


PROCESS PART 

DEFINE TTF AS A REAL VARIABLE 
DEFINE RT AS A REAL VARIABLE 
DEFINE DT AS A REAL VARIABLE 
DEFINE RC AS A REAL VARIABLE 
DEFINE DC AS A REAL VARIABLE 
DEFINE TC AS A REAL VARIABLE 
UNTIL TIME. V >= RUN. LENGTH 
DO 


''TIME TO FAILURE 
''REPAIR TIME 
''DELAY TIME 

''REPAIR COST PER REPLACEMENT 
"DELAY COST PER INCIDENT 
''TOTAL COST PER INCIDENT 


ADD 1 TO NOW.WORK 

WORK WEIBULL. F(FAIL. SHAPE,FAIL. SCALE,SEED1). MILES TIME TO FAIL 
LET TTF = WEIBULL. F(FAIL. SHAPE,FAIL. SCALE,SEED1) " 

IF TTF > 15450 

'' PRINT 1 LINE WITH NO. FAIL THUS 
'' ***** skeD REPLACE EARLY 

TTF = 15450 

WORK TTF .MILES "UNTIL FAILURE TIME 
SUBTRACT 1 FROM NOW. WORK 
" PRINT 1 LINE WITH NO. FAIL THUS 
'' ***** part REPLACED EARLY 

ADD 1 TO REPL.NO 

LET RT = LOG. NORMAL. F(REP. MU,REP.SIG.SEED3) 

LET RC = PART. COST + (MN. HR. COST * RT) 
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ADD RC TO REP. COST 


WAIT RT HOURS ''TIME TANK UNAVAILABLE 
LET DC = (STAND-BY. COST/24) * RT 
ADD DC TO PEN. COST 
LET TC = RC + DC 
ADD TC TO TOT. COST 
CYCLE 
ALWAYS 

WORK TTF .MILES "UNTIL FAILURE TIME 
SUBTRACT 1 FROM NOW. WORK 
ADD 1 TO NO. FAIL 

LET RT = LOG. NORMAL.F(REP.MU,REP.SIG.SEED3) ’’ 

LET RC = PART. COST + (MIL HR. COST * RT) 

ADD RC TO REP. COST 

" WAIT LOG. NORMAL. F(DELAY. MU,DELAY. SIG.SF.ED2) DAYS DELAY TIME 
LET DT = LOG.NORMAL. F(DELAY. MU,DELAY. SIG.SEED2) 

WAIT DT DAYS "TIME TANK UNAVAILABLE 
LET DC = STAND. BY. COST * DT 
ADD DC TO PEN.COST 
LET TC = RC + DC 
ADD TC TO TOT.COST 
LOOP 

END 


PROCESS STOP.SIM 

LET T. FAIL. .= NO. FAIL 
LET PRCNT. UP. . = AVAIL 
LET R. COST. . = REP. COST 
LET D. COST. . = PEN. COST 
LET T. COST. . = TCT. COST 

"PRINT 1 LINE WITH N AND TIME. V THUS 
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1 * * **_ ** 

''PRINT 1 LINE WITH NO. FAIL THUS 
'* ***** FAILED 

’’PRINT 1 LINE WITH AVAIL*100/N THUS 
" ***.**** PERCENT AVAIL 

''PRINT 1 LINE WITH REP. COST THUS 
'* ***.**** COST REP 

''PRINT 1 LINE WITH PEN. COST THUS 
•» ***, **** COST PEN 

''PRINT 1 LINE WITH TOT. COST THUS 
'' ***. **** COST TOTAL 

''******** 

” PRINT 1 LINE WITH N,MILES. N.DOWN. N , NO. FAIL AND AVG. DOWN. T THUS 
' ' * **, ** ** ' ** ** ** ** 

TIME. V = 0.0 
'' DOWN. N = 0. 0 
' ' MILES. N = 0. 0 
NO.FAIL = 0 
REP. COST = 0.0 
PEN. COST = 0.0 
TOT. COST = 0.0 

DESTROY THIS PART 
' PRINT 1 LINE WITH TIME.V THUS 
' ** ** DAYS 

' PRINT 1 LINE WITH DOWN. T THUS 
' **.** DAYS DOWN 

' PRINT 1 LINE WITH MILES. T THUS 
’ **,** MILES 

’ PRINT 1 LINE WITH NO.FAIL THUS 
' ***** FAILED 

IF N >= 5000 
CALL LAST. RUN 
ALWAYS 
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END 

ROUTINE LAST. RUN 

PRINT 1 LINE WITH N THUS 
FOP. **** RUNS 

PRINT 1 LINE WITH REPL. NO THUS 
**** PARTS REPLACED EARLY UNDER POLICY. 

PRINT 1 LINE WITH AVG.T. FAIL AND SIG. T. FAIL THUS 

AND AVERAGED ******* FAILURES WITH STD. DEV= ****.** IN YEARS. 

PRINT 1 LINE WITH AVG. R. COST AND SIG. R. COST THUS 

AND AVERAGED *******.** REPAIR COST WITH STD. DEV= ******.** 

PRINT 1 LINE WITH AVG. D. COST AND SIG.D.COST THUS 

AND AVERAGED *******. ** PENALTY COST WITH STD. DEV= ****** ** 

PR I ITT 1 LINE WITH AVG. T. COST AND SIG. T. COST THUS 

AND AVERAGED *******. ** TOTAL COST WITH STD. DEV= ******. ** 

PRINT 1 LINE WITH PRCNT.UP*100/N THUS 
***.**** PERCENT AVAILABLE 
STOP 
END 
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APPENDIX I. SYSTEM SIMULATION 


Code: SIMSCRIPT 11.5 
Programmer: J. Wilhelm 
Date: July 1990 

’' TANK RAM MODEL 

PREAMBLE 

NORMALLY, MODE IS UNDEFINED 
DEFINE . MILES TO MEAN HOURS 
PROCESSES INCLUDE TANK. DOWN AND STOP.SIM 
EVERY PART 

BELONGS TO A PART. SET 
AND HAS A FAIL. SHAPE 
AND HAS A FAIL.SCALE 
AND HAS A DELAY. MU 
AND HAS A DELAY. SIG 
AND HAS A SEED1 
AND HAS A SEED2 

DEFINE FAIL.SHAPE AS A REAL VARIABLE 
DEFINE FAIL. SCALE AS A REAL VARIABLE 
DEFINE DELAY. MU AS A REAL VARIABLE 
DEFINE DELAY.SIG AS A REAL VARIABLE 
DEFINE SEED1 AND SEED2 AS INTEGER VARIABLES 

THE SYSTEM OWNS THE PART.SET 

DEFINE NO.OF.PARTS AND NO. FAIL AS INTEGER VARIABLES 
DEFINE N AS A INTEGER VARIABLE 
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DEFINE T. FAIL AS A INTEGER VARIABLE 
DEFINE RUN. LENGTH AS A REAL VARIABLE 
DEFINE MILES. T AS A REAL VARIABLE 
DEFINE MILES. N AS A REAL VARIABLE 
DEFINE DOWN. T AS A REAL VARIABLE 
DEFINE DOWN.N AS A REAL VARIABLE 
TALLY AVG. MILES. T AS THE MEAN, 

SIG. MILES. T AS THE STD. DEV OF MILES. T 
TALLY AVG. DOWN. T AS THE MEAN, 

SIG. DOWN. T AS THE STD. DEV OF DOWN.T 
TALLY AVG. T. FAIL AS THE MEAN, 

SIG. T. FAIL AS THE STD. DEV OF T. FAIL 

END 

MAIN 

LET HOURS. V = 6.5373 "MILES 
CALL READ. DATA 
FOR N = 1 TO 5,DO 
CALL INITIALIZE 

ACTIVATE A STOP. SIM IN RUN. LENGTH DAYS 

START SIMULATION 

LOOP 

END 

ROUTINE INITIALIZE 
"DEFINE I AS A INTEGER VARIABLE 
CREATE A TANK. DOWN 

' ' FOR I = 1 TO NO. OF. PARTS 
' ' DO 

'' STARTER 
CREATE A PART 

LET FAIL.SHAPE(PART) = 1.1999 
LET FAIL.SCALE(PART) = 4689.7 
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LET DELAY. MU (PART) = .59102 

LET DELAY. SIG(PART) = . 79943 

LET SEED1(PART) = 1 

LET SEED2(PART) = 2 

FILE THIS PART IN THE PART. SET 

ACTIVATE THIS PART NOW 

’' FUEL NOZZLE 

CREATE A PART 

LET FAIL. SHAPE (PART) = 1.3384 

LET FAIL. SCALE(PART) = 4438.7 

LET DELAY. MU (PART) = .69437 

LET DELAY. SIG(PART) = .67156 

LET SEED1(PART) = 3 

LET SEED2CPART) = 4 

FILE THIS PART IN THE PART. SET 

ACTIVATE THIS PART NOW 

'' DIST BOX 

CREATE A PART 

LET FAIL. SHAPE(PART) = 1.5315 

LET FAIL.SCALE(PART) = 6858.6 

LET DELAY. MU (PART) = .57369 

LET DELAY. SIG(PART) = .93421 

LET SEED1(PART) = 5 

LET SEED2(PART) = 6 

FILE THIS PART IN THE PART. SET 

ACTIVATE THIS PART NOW 

'' TRANS 

CREATE A PART 

LET FAIL. SHAPE(PART) = 1.4455 
LET FAIL. SCALE(PART) = 7479.5 
LET DELAY.MU (PART) = .15064 
LET DELAY. SIG(PART) = .94184 
LET SEEDl(PART) = 7 
LET SF.ED2(PART) = 8 
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FILE THIS PART IN THE PART.SET 
ACTIVATE THIS PART NOW 
'' LINK 
CREATE A PART 

LET FAIL. SHAPE (PART) = 1.2801 

LET FAIL. SCALE(PART) = 10239.0 

LET DELAY. MU (PART) = .3510 

LET DELAY. SIG(PART) = 1.1755 

LET SEED1(PART) = 9 

LET SEED2(PART) = 1 

FILE THIS PART IN THE PART. SET 

ACTIVATE THIS PART NOW 

'' GRIP 

CREATE A PART 

LET FAIL.SHAPE(PART) = 1.5539 
LET FAIL. SCALE(PART) = 4369.9 
LET DELAY.MU (PART) = .60291 
LET DELAY. SIG(PART) = .87986 
LET SEED1(PART) = 2 
LET S£ED2(PART) = 3 
FILE THIS PART IN THE PART.SET 
ACTIVATE THIS PART NOW 
LET FAIL. SHAPE(PART) = 1.0 
LET FAIL. SCALE(PART) = 5.0 
LET DELAY. MU (PART) =2.4 
LET DELAY. SIG(PART) = 1.0 
LET SEEDl(PART) = 1 
LET SEED2(PART) = 1 

LOOP 

ACTIVATE A TANK. DOWN NOW 

ACTIVATE A UP. TANK NOW 

LET FAIL.SCALE(1) = 5. 0 


END 
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ROUTINE READ. DATA 


” PRINT 1 LINE THUS 
'' HOW MANY TANKS? 

'' READ NO.OF.TANKS 
" NO.OF. TANKS = 1 

" PRINT 1 LINE THUS 
'' HOW MANY PARTS PER TANK? 

'' READ NO. OF. PARTS 

NO.OF.PARTS = 6 
'' PRINT 1 LINE THUS 
" HOW LONG SHALL WE RUN (IN DAYS)? 
'' READ RUN.LENGTH 

RUN. LENGTH = 7300.0 

' ' DOWN. T = 0. 0 
! ' MILES. N = 0. 0 
' ' LET FAIL. SCALE. .=5.0 
'' FAIL. SCALE(2) = 5. 0 
END 


PROCESS PART 

UNTIL TIME. V >= RUN. LENGTH 
DO 

WORK WEIBULL. F(FAIL. SHAPE.FAIL. SCALE,SEED1). MILES 
REACTIVATE THE TANK. DOWN NOW 

WAIT LOG.NORMAL.F(DELAY. MU,DELAY. SIG.SEED2) DAYS 
REACTIVATE THE TANK.DOWN NOW 
LOOP 

'' IF TIME.V <= RUN. LENGTH 

'' REMOVE THIS PART FROM THE PART.SET 

" ALWAYS 

" SUSPEND 

END 


"TIME TO FAIL 
' REPAIR TIME 
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PROCESS TANK. DOWN 


DEFINE FAIL. T AS A REAL VARIABLE 
FAIL. T - 0. 0 

DEFINE REPAIR. T AS A REAL VARIABLE 
REPAIR.T = 0.0 
UNTIL TIME. V >= RUN. LENGTH 
DO 

SUSPEND "AWAITING A SUBASSEBLY FAILURE 
LET FAIL. T = TIME. V 
'' INTERRUPT UP. TANK 

FOR EACH PART IN THE PART. SET, 

WITH STA. A(PART) = 1 "i.E. ,IT IS OPERATING 
DO 

INTERRUPT THIS PART 
LOOP 

ADD 1 TO NO. FAIL 

SUSPEND "AWAITING REPLACEMENT OF SUASSEMBLY 
LET REPAIR. T = TIME. V - FAIL. T 
LET DOWN.N = DOWN. N + REPAIR. T 
LET DOWN. T = DOWN. N 

" LET MILES. T. ,= 1 * (TIME. V - DOWN. T) 

'' RESUME UP. TANK 

FOR EACH PART IN THE PART. SET, 

WITH M. EV. S( PART) <> 1 "I.E. ,IT IS NOT SCHEDULED 
DO 

RESUME THIS PART 
LOOP 
LOOP 

SUSPEND "AWAITING END OF SIMULATION 


END 


no 






PROCESS STOP.SIM 


LET MILES. N. .= 6.5373 * (TIME. V - DOWN. T) 

LET MILES. T. .= 6.5373 * (TIME.V - DOW. T) 

LET T. FAIL.. = NO.FAIL 
'' PRINT 1 LINE WITH N AND TIME. V THUS 

* 1 A A*. A A 

" PRINT 1 LINE WITH N,MILES. N.DOWN. N , NO. FAIL AND AVG. DOWN. T THUS 

' ' A AA,A* A*.AA AA AA, AA 

TIME. V = 0. 0 
DOWN. N = 0. 0 
MILES. N = 0. 0 
NO. FAIL = 0 

IF M. EV. S(TANK. DOWN) = 1 

REMOVE TANK.DOWN FROM EV. S(I. TANK.DOWN) 

'' DESTROY TANK. DOWN 
ALWAYS 

FOR EACH PART IN THE PART. SET, 

WITH M.EV. S(PART) = 1 
DO 

REMOVE THIS PART FROM EV. S(I.PART) 

LOOP 

FOR EACH PART IN THE PART.SET 
DO 

REMOVE THIS PART FROM THE PART. SET 
DESTROY THIS PART 
LOOP 

" REMOVE PART FROM EV.S(I.PART) 

" PRINT 1 LINE WITH N.EV.S (I.PART) THUS 
' ' * NUM ON EV SET 
DESTROY TANK. DOWN 
' ' DESTROY PART 

” PRINT 1 LINE WITH TIME.V THUS 
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'’ **, ** DAYS 

'' PRINT 1 LINE WITH DOWN. T THUS 
'' **.** DAYS DOWN 

'' PRINT 1 LINE WITH MILES. T THUS 
'' **.** MILES 

'' PRINT 1 LINE WITH NO. FAIL THUS 
'' ***** FAILED 

IF N >= 5 

CALL LAST. RUN 
''ELSE 
'' RETURN 
ALWAYS 
END 

ROUTINE LAST. RUN 

PRINT 1 LINE WITH N,AVG. MILES.T AND SIG.MILES. T THUS 
FOR **** RUNS TANKS AVERAGED ****** WITH STD. DEV= **. ** (MILES) 
PRINT 1 LINE WITH AVG. DOWN. T AND SIG. DOWN. T THUS 
THEY WERE DOWN AN AVERAGE OF **. ** WITH STD. DEV= **.** (DAYS) 
PRINT 1 LINE WITH AVG. T. FAIL AND SIG. T. FAIL THUS 
AND AVERAGED **.** FAILURES WITH STD. DEV= **.** IN TWO YEARS. 
STOP 
END 


V, 


< 


( 


L 
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